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Abstract 

When nonlinear dynamical systems are coupled, depending on the intrinsic 
dynamics and the manner in which the coupling is organized, a host of novel 
phenomena can arise. In this context, an important emergent phenomenon 
is the complete suppression of oscillations, formally termed amplitude death 
(AD). Oscillations of the entire system cease as a consequence of the inter- 
action, leading to stationary behavior. The fixed points that the coupling 
stabilizes can be the otherwise unstable fixed points of the uncoupled system 
or can correspond to novel stationary points. Such behaviour is of relevance 
in areas ranging from laser physics to the dynamics of biological systems. In 
this review we discuss the characteristics of the different coupling strategies 
and scenarios that lead to AD in a variety of different situations, and draw 
attention to several open issues and challenging problems for further study. 

Keywords: Amplitude Quenching, Interaction, Synchronization, 
Fixed-point solution, Control, Network, Bifurcation 



Email addresses: gsaxenal@physics . du. ac . in (Garima Saxena), 
awadhesh@physics.du.ac. in (Awadhesh Prasad*), rrOuohyd . ernet . in (Ram 
Ramaswamy) 

* Corresponding Author: Fax +11-2766 7061; Phone +11-2766 7725/Extn. 1354. 



Preprint submitted to Physics Reports 



October 1, 2012 



Contents 



1 Introduction 

2 Scenarios for Amplitude Death 6 

2.1 Mismatched oscillators [7 

2.2 Delay interaction 

2.3 Conjugate coupling 

2.4 Dynamic coupling 

2.5 Nonlinear Coupling 

2.6 Linear augmentation, and other strategies 

2.7 Velocity Coupling 

2.8 Routes to AD 

2.8.1 Hopf bifurcation 

2.8.2 Saddle-node Bifurcation 

2.8.3 Direct transition 

3 Targeting and Control 

3.1 Steady-states through nonlinear coupling 

3.1.1 Linear Augmentation 

3.2 Reprieve: Avoiding fixed points 

3.2.1 Gradient coupling 

4 Characterization 

4.1 The phase-flip 

5 Amplitude death on networks 

6 Experiments and applications 

7 Summary 

8 Acknowledgments 

9 References 



2 



1. Introduction 

The nature of the dynamics in coupled nonlinear systems is a central 
issue of interest in a number of areas of scientific enquiry. A variety of 
natural phenomena-ranging from the motion of pendulums or springs, to 
lasers, fluid flow, ecological systems, populations, and so on-are modeled as 
oscillators. These may be linear or nonlinear, depending on the context, 
and as is well known, the motion can be periodic or chaotic. When systems 
that are in isolation capable of complex behaviour are coupled, a host of 
novel phenomena can be seen. These depend both upon the properties of 
the isolated systems — the nature of the nonlinearity, whether the motion is 
chaotic or regular, what are the equilibria, and so on — as well as the manner 
in which their coupling is organized. 

An important emergent phenomenon in this context is the suppression 
of oscillations, formally termed amplitude death (AD): as a consequence of 
the interaction oscillations of the entire system cease, leading to stationarity 
0, 0, El- Such behaviour is emergent in the sense that the isolated or uncou- 
pled systems do not exhibit stationary dynamics. In a formal sense, ampli- 
tude death is an instance of a more general phenomenon that includes both 
the cessation of oscillations as a consequence of coupling, as well the suppres- 
sion of amplitude variations, namely the conversion of chaotic oscillations to 
periodic or quasiperiodic dynamics. In a sense, this is the stabilization or 
creation of "simpler" attractors through the coupling. 

When oscillations stop, there are two possibilities for AD. If the cou- 
pled systems have exactly one equilibrium, then the occurrence of amplitude 
death implies that this equilibrium becomes asymptotically stable. However, 
the coupled systems can have more than one stationary state and then one 
can have the stabilization of new fixed points, namely those that are not 
stable (or may not even exist) in the uncoupled system. 

A general setting in which we discuss AD in this review is a configuration 
of N coupled nonlinear oscillators that is specified by the equation 

X l = F 4 (X 4 )+eG 4 (X) z = l,...,iV. (1) 

where 

1 N 

G*(X) = — J^-H^X^r), * = (2) 
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Here X, is the 7?v dimensional vector of dynamical variables for the ith oscil- 
lator, and the Fj's specify their evolution equations. The oscillators may be 
identical or distinct, and are coupled to other oscillators, as specified via the 
function Gj. Ki is a normalization factor, e is the coupling strength, and the 
connection topology is specified by the connectivity matrix, whose elements 
Aij are 1 or 0, depending on whether oscillators i and j are coupled with 
each other or not. The actual coupling is specified by the term H(Xj, Xj, r) 
which is a function of Xj(t) and Xj(£ — r). The time-delay r can be fixed, 
or can have a more complex dependence on time. Most of the major results 
discussed in this review are illustrated with two coupled oscillators, N = 2. 

Shown in Fig. [I] is a schematic of the basic phenomenology. Consider two 
oscillators coupled by some function G. When the coupling strength e = 
the motion of individual subsystems are oscillations (whether chaotic (C), 
periodic (P) or quasiperiodic (QP)); a simple periodic case is illustrated in 
Fig. [T^a). When the coupling is switched on, above a specific strength when 
amplitude death occurs, the transient dynamics is essentially that of damped 
oscillations, as shown in Figs. [TJ (b) and (c). As can be seen, although this 
is an amplitude phenomenon, there could be additional effects that involve 
the phases of the oscillators, and indeed, any analytic understanding of this 
phenomenon must include a discussion of both the amplitudes and the phases 
of the interacting systems. 

There are some terminological issues that need to be clarified. In this ar- 
ticle we will term the cessation of oscillations, regardless of what fixed points 
are stabilized asymptotically, as amplitude death, even though the asymp- 
totic steady states in the different subsystems may not be the same (as for 
example in Fig. He)) and can be different from the null state, X* = 0. This 
heterogeneous case has been termed "oscillation death" [l], 0, [i[ or the Bar- 
Eli Effect jsl), keeping the term amplitude death for the homogeneous case, 
when a simple transformation can shift the fixed point to the origin 0, 0] (as 
for example in Fig. QJb)). What these fixed points are will naturally depend 
on the internal structure of the individual oscillators, whether they are iden- 
tical or mismatched, whether the coupling is symmetric or asymmetric, and 
indeed what the nature of the coupling is, namely whether the interactions 
are diffusive or nonlinear. 

Following the initial work by Bar-Eli [l[ AD has been studied quite exten- 
sively jsl, Is), H, 10 1 . Different situations where amplitude death occurs have 



been described, and it is now known that the phenomenon can be seen when 
the interacting systems are identical @, EH 12, 13], mismatched [3, 0], when 
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Figure 1: (Color online) Schematic of the dynamics for a typical case of two oscillators 
with a given coupling function G, showing (a) the phase space dynamics of the uncoupled 
oscillators, namely when e = 0, and the transient motions of the two subsystems (solid 
and dashed lines for respective oscillators) as a function of time, leading to either (b) 
the same (homogeneous) or (c) the distinct (heterogenous) steady states respectively for 
appropriate coupling strength, e 7^0. 
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the coupling is dynamical 14 1. or nonlinear (nondiffusive) [15[. We review 
these various situations and also discuss the general mechanism for AD in 
Sec. m 

)ractical situations AD, is desirable, as for example in laser 
16|, [l7|, [l8|, yjj, |2(J, |2l| where a constant output is needed and 
fluctuations should be suppressed. There are also other situations where 



In many 
applications 



oscillations need to be maintained, as for instance in brain functioning [22 



23 



These different requirements suggest that control strategies to either to 
achieve or to avoid AD in coupled systems are necessary, and in Sec. [3] we 
review some control methods that have been employed to either target or 
avoid steady states in coupled systems. 

As indicated above, understanding the transient dynamics is important 
in describing AD since this state is asymptotically featureless. Transient dy- 
namics can be significant in many applications that are restricted to finite 
times, as for example in ecology 24J. Furthermore, since the different sub- 
systems are coupled, a study of the synchronization properties 25] is also 
of interest. In a number of examples it has been noted that the individual 
systems, while being synchronized, can nevertheless be either in-phase or 
out of phase, with a transition between these states at a specific value of the 



coupling [26|, |27J, |28| . We review this phase-flip transition and its associated 
behaviour in Sec. [U 

The occurrence of AD in networks of coupled oscillators 



15, 29, 30, 31 



32|, [33|, [34j, [35[ is another topic of considerable current relevance. Studies have 



been carried out for a variety of network topologies, ranging from systems 
coupled on a ring 
scale-free networks 



36, 37 



381 and 



as well as in the so-called small world 
These are discussed in Sec. [5] 
Our principal aim in this review on amplitude death is to stress the 
importance of this phenomenon and its applications in real situations. The 
major experimental results and possible applications are discussed in Sec. [6j 
Finally, the review concludes in Sec. [7] with a summary. 



2. Scenarios for Amplitude Death 

To begin with, we consider the several known scenarios within which AD 
is known to occur, as well as constraints and necessary conditions. One aspect 
of interest is in identifying the different types of interactions that facilitate 
AD. This is of particular relevance in natural systems where the parameter 
of individual systems may not be accessible, and amplitude death can be 
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achieved mainly by selecting appropriate forms of the interaction. It is also 
possible that for a given form of the interaction, AD occurs when internal 
parameters are varied. As a parallel, consider a related phenomenon, the 
synchronization of two systems. There are specific forms of interaction that 
cause systems to synchronize, and for a specific form of coupling — diffusive, 
say — synchrony is possible only when the coupling strength is sufficiently 
strong. 

The Landau-Stuart limit cycle oscillator system 0, Ho|, 41 1 has been an 



ideal model for the exploration of several of the scenarios we discuss here. 
This is a two-dimensional oscillator given by the evolution equation 

Z = [1 +itu - \Z\ 2 ]Z, (3) 

where the complex variable Z = x + iy. As can be seen, the origin is an 
unstable fixed point, and the unit circle is an attractor for positive frequency, 
oo. The relative simplicity of the evolution equations makes the analysis 
tractable, including the case of time-delayed interactions. 

In this section the case of two oscillators is considered, and some numerical 
results are presented first. The stability of the fixed points can be estimated 
by computing the eigenvalues of the Jacobian matrix. When all eigenvalues 
have negative real parts, the entire spectrum of Lyapunov exponents 42| is 



negative: this is the state of AD. In the numerical results presented here the 



standard Runge-Kutta scheme [43j was used. 



2.1. Mismatched oscillators 

An early observation of amplitude quenching was reported by Crow- 



ley and Field [44J who experimentally coupled two chemical oscillators, the 
Belousov-Zhabotinski system electrically. The coupling essentially involved 
the mass transfer of a single species, the Ce 3+ ions, and caused the composite 
system to reach a state wherein one of the systems was quenched to a steady 
state, while the other maintained oscillations. This experiment served as 
the starting point for the study of this behaviour by Bar-Eli [l[ who showed 
that two interacting model continuous stirred tank reactors can stop oscil- 
lations and arrive at steady states if coupled diffusively. Consider the cell 
Brusselator model [if, defined as 

x = -(B + l)x + x 2 y + A 

y = Bx — x 2 y (4) 
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where x and y are the dimensionless concentrations of two chemical species 
[il, [l[ while the (positive) parameters A and B are the rate constants for the 
production of x and y. This system has one steady state, x* = A, y* = B/A 
which is unstable if A 2 < (B — 1). The oscillating behavior of this system is 
shown in Fig. H](a) for parameters A = 5 and B = 35. When two such cells 
are coupled diffusively |l| then the dynamical equations take the form (the 
variables of the two systems are distinguished by subscripts) 

x\ = -(Bi + l)xi + x\yi + A 1 + e(x 2 - xi) 

y\ = B x x x - x\y x + e(y 2 - yi) (5) 

with corresponding equations for subsystem 2. The coupling terms e(x2 — x x ) 
and e{y2—y\) arise from material transfer between the cells through diffusion, 
and the coupling strength e corresponds to the rate of this transfer. The 
largest two Lyapunov exponents at A x = A2 = 5, B x = B2 = 35 as a function 
of the coupling strength in Fig. [2](b). The transition to AD can be identified 
as the point when all Lyapunov exponents become nonpositive (and which 
then correspond to the real parts of the eigenvalues of the fixed point-see 
Sec. H]). Both the oscillators then settle to different steady states. 

Transient trajectories for both oscillators are shown in Figs. [2](c) and (d), 
demonstrating that this is a case of oscillation death: the two fixed points 
(xl = 0.34, y* = 22.36) and (x* 2 = 9.66,^ = 3.72) are distinct, and different 
from the fixed point of the uncoupled systems, namely x* = 5, y* = 7. 

Similar behavior was found to exist in other well-known chemical oscil- 
lators [3, Hi|, including the Noyes-Field-Thompson model, the Oregonator, 
and the Field-koros-Noyes models for the Belousov-Zhabotinski reaction, 



the model of first-order decomposition autocatalysis [46J and the Lotka- 
Volterra model for prey-predator interactions. These systems differ signif- 
icantly showing that the occurrence of amplitude death is fairly general. 
This extensive numerical analysis [l| triggered a serious exploration of the 
phenomenon of amplitude death from both a theoretical and experimental 
point of view. 

Analysis by Aronson, Ermentrout, and Kopell [s| provided a deeper math- 
ematical understanding of this phenomenon. They considered the following 
system of two Landau-Stuart oscillators, 

Z x = [l + iWi-lZtflZi + etZa-Zx], 

Z 2 = [1 + iu 2 - \Z 2 \ 2 ]Z 2 + e[Z x - Z 2 ] (6) 
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Figure 2: (Color online) (a) The trajectory of a single Brusselator, Eq. (j4]), at parameters 
A = 5 and B = 35 in phase space, (b) The largest two Lyapunov exponents as a function 
of coupling strength e for two coupled identical Brussclators, Eq. at parameters A\ = 
A2 = 5 and Bi = B2 — 35. (c) And (d) show the transient trajectories of both oscillators 
as function of time at coupling strength e = 0.5. 
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Figure 3: The phase diagram in the parameters: frequency mismatch (Aw) and coupling 
strength (e) of Eq. ((6|). The shaded and unshaded regions correspond to the periodic (P) 
and amplitude death (AD) motions respectively. For detail see Ref. [f|. 



where the notation is obvious. When uncoupled, both oscillators have an 
unstable fixed point at Z\ 2 = 0, and the unit circles |Z li2 1 = 1 are limit 
cycles. Depending on the mismatch Au = \oji — a>2|, this system shows 
different types of dynamics j^J. Shown in Fig. [3] is a schematic phase- 
diagram, indicating the different dynamical states: the shaded and blank 
regions correspond to periodic motion and AD respectively. For a specific 
value of the mismatch, Au = 4, the largest two Lyapunov exponents are 
shown in Fig. HI as a function of the coupling. The transient dynamics in 
the AD region is shown in the inset; note that this is a case of homogenous 
steady states and both oscillators settle onto the (now stabilized) fixed point, 
the origin Z* = 0. Linear stability analysis in the neighborhood of the origin 
shows that amplitude death occurs for e > 1 and Au > 2y/2e — 1. 

The sign of u determines the sense of rotation in phase. In Eq. fl6]) the fre- 
quencies are both positive so that in the phase plane the direction of motion 
for both oscillators is the same. In this case, the systems are termed co- 
rotating, while when the frequencies differ in sign, they are counter-rotating 



47 1 . This latter situation is one of parameter mismatch, and in recent works 



471 1481 . 149] . the characteristics of AD in such cases have been discussed; see 
also Sec. HI 

Parameter mismatch is thus one of the major causes of amplitude death, 
and has been extensively studied both analytically and experimentally [H, B 
S S |47 . 
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Figure 4: (Color online) The largest two Lyapunov exponents as a function of coupling 
strength, at Aoj = 4 {w\ = 4 and cj 2 =8), Eq. (|6]). Inset figure shows in-phase synchronous 
behaviour of transient trajectories of both the oscillators as a function of time at coupling 
strength e = 1.5 , x\ (solid-black line) and Xi (dashed-red line). 



2.2. Delay interaction 

In the above systems, the interaction is considered to be instantaneous, 
namely the coupling terms involve the variables of the two subsystems at the 
same time, t. In many physical systems, the coupling involves time-delay, 
and this was considered by Reddy, Sen and Johnston [7J who investigated the 
collective dynamical behavior of limit-cycle oscillators interacting diffusively 
via time-delayed coupling. Here AD was achieved even in the absence of 
mismatch, and was demonstrated in an experiment involvin g el ectronic cir- 
cuits. As was shown through both analysis and numerics fK l50l. 51. 52. 53 



there is a significant region in the parameter space of coupling strength and 
time delay where AD occurs. 

The delay-coupled Landau-Stuart system has the evolution equations 

Z 2 = [l + iu2- \Z^\Z 2 + s[Z x {t-T) -Z 2 ], (7) 

where the output of one system takes a time r to influence the other, as 
might occur in spatially separated systems when there is a finite speed of 
signal transmission. 

Consider identical systems, u\ = uo 2 = 6. Shown in Fig. \5\is the largest 
two Lyapunov exponents as a function now of the time-delay r for fixed 
coupling strength e — 2. This clearly shows the regime of AD when all Lya- 
punov exponents are negative (see details in Sec. 12 . 8 j) . Transient trajectories 
(shown in the inset) for both oscillators indicates the cessation of oscillations 
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Figure 5: (Color online) The largest two Lyapunov exponents as a function of coupling 
strength for coupled identical oscillators for Aw = (w\ = 6 and o>2 = 6), (Eq. ([?])). 
Inset figure: completely synchronized transient trajectories, x\ (solid-black line) and X2 
(dashed-red line), as a function of time at time-delay r = 0.255 @, Upl- 



and the stabilization of the homogeneous steady state, Z\ 2 = 0. The AD 
in such delayed coupled systems has a resonance-like character, occurring in 
repeated "islands" in the parameter space of e and r 0, [Hq|, 51, 52, 53 . 

This scenario for the occurrence of AD in delay-coupled systems is quite 
general and does not depend on the nature of the dynamics in the uncoupled 



systems. As was shown by Prasad [26|, even chaotic oscillators can be sta- 
bilized in this manner. The importance of time delay interactions has been 



highlighted in [54 



Given a finite transmission speed for signals from one subsystem to reach 
the other, a natural extension is to consider the situation when the time- 
delay itself varies. This is particularly important when there are stochastic 
effects that are being modeled by the delay itself. Atay [l3| considered the 
effect of a distribution of time-delays and showed that even a small spread 
in the delay distribution can greatly enlarge the set of parameters for which 
amplitude death occurs. 

The expansion of the AD region with a distribution of delays is a impor- 
tant effect with a wide range of applications 13j, particularly since if the 



variance of the distribution is above a threshold, this AD region can become 
unbounded. Again, the effect is quite general, occurring for arbitrary spread 
in the delay, for different forms of the distribution, and for arbitrary numbers 



12 



of oscillators [13|, |55| . Closely related to the case of distributed delays is the 



situation where systems respond to cumulative signals, namely by integrating 
information received over an interval in time. This occurs when systems have 
a finite intrinsic response time and also causes the region of AD to extend 
indefinitely 



56. 57 



Since time-delay interactions are common in experiments, various aspects 
of both fixed and distributed delays, including a deterministic time- varying 
have been studied. One application is to the stabilization 
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interaction 



of fixed points 



59 
60 



61, 62 



Velocity delayed coupling has been also explored 
in the context of amplitude death [63j, |64| (see Sec. 12. Tj) . Mixed interaction 
schema, with instantaneous interactions in one system and time-delay in the 



other 65[ or unequal time-delays in different directions (to model spatial 



66|, [ 67j and one-way ring time-delay |37J have been studied in 



heterogeneity 
recent works. 



2.3. Conjugate coupling 

A novel context for amplitude death is when systems are coupled via 
so-called conjugate or dissimilar variables. With such coupling, the neces- 
sity for systems being either mismatched or having time-delayed interactions 
can be dropped, and AD occurs in identical systems with instantaneous cou- 
pling. This was first studied in a model system by Karnatak et al. 11| who 
considered two Landau-Stuart oscillators with the evolution equations 



± 1 = P lXl - utyt + s(y 2 - xi) 
2/i = PiVi + + £(x 2 - ?/i ) 



(8) 



where Pi = 1 — \Zi\ 2 , and similar equations can be written for subsystem 2. 
Note that here the y variable of the second oscillator is diffusively coupled 
to the x variable of the first oscillator and vice versa. 

Coupling via conjugate variables is natural in a variety of experimental 
situations, and indeed the above system was inspired by the experiments 
of Kim and Roy 16|, llTJ on coupled semiconductor laser systems where the 
photon intensity fluctuations from one laser were used to modulate the in- 
jection current of a second identical laser and vice versa. Indeed it would 
not be possible to couple two lasers through their intensity fluctuations di- 
rectly. As it happens, in a model of the semiconductor laser, the current and 
the fluctuations are, in a sense, conjugate variables that have constitute the 
chaotic oscillator system 16j, |17|. The largest three Lyapunov exponents for 
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Figure 6: (Color online) The largest three Lyapunov exponents as a function of coupling 
strength at Auj = (u>i = 6 and ui2 = 6), Eq. ((5J. Inset figure: the transient trajectories, 
X\ (solid-black line) and xi (dashed-red line), as a function of time at coupling strength 
e = 50. 



the above system is shown in Fig. |6l and as can be seen, there is a wide range 
of coupling strength where all Lyapunov exponents are negative, indicating 
stabilization of the fixed point. A transient trajectories are shown in the 
inset. 

A detailed analytical condition for the occurrence of amplitude death can 



be found in [68|. With conjugate coupling, AD can be anticipated by rec- 



ognizing that there is a strong analogy between time-delayed and conju- 



gate variables: as a matter of fact, the Takens' embedding theorem [69 
asserts that the topological properties of a dynamical system that is recon- 
structed through delay variables match those of the true system for appro- 
priate choices of embedding dimension and time delay. Thus using conjugate 
variables is akin to using time-delay coupling, and this gives rise to regimes 
of AD. 

In addition, there are other new features such as the suppression of chaos, 
and riddling 0, Q [sj due to conjugate coupling. In addition, this is 



particularly suited to experiments and both in electrochemical systems [71 



as well as in electronic circuits [72[|, this coupling can be realized. 
2.4- Dynamic coupling 



Konishi [14J has proposed another type of interaction which, in contrast to 
the cases discussed above has an evolving or "dynamic" coupling. Consider 
the Landau-Stuart system 

±t = PtXi - uiyi + e(ui - x{) 
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Figure 7: (Color online) The largest two Lyapunov exponents as a function of coupling 
strength, at Au; = (u>i = 6 and ui2 = 6), Eq. (|9]). Inset figure: transient trajectories, 
X\ (solid-black line) and xi (dashed-red line), as a function of time at coupling strength 
£ = 3 0. 



yi 



PiUi + UJ1X1 



(9) 



with similar equations for subsystem 2, and where Pi = 1 — \Zi\ 2 as usual. 
The variables Ui follow linear dynamics in absence of the x/s so that the 
coupling has its own nontrivial dynamics. Shown in Fig. [7] is the largest two 
Lyapunov exponents; the occurrence of amplitude death is evident, and the 
transient dynamics are shown in the inset. Note that this instance of ampli- 
tude death is in identical systems. In addition to estimating necessary and 
sufficient conditions for amplitude death to occur, this form of the dynamics 



has lent itself to extensive analytic study [73|, |74|, [75 



2.5. Nonlinear Coupling 

In the preceding examples, the form of the coupling has been taken to 
be linear or diffusive. An interesting extension that uses nonlinear coupling 
makes it possible to not only achieve AD, but also to specify the steady state 
of the coupled system 15|, |2l|. We will illustrate the targeting of specific 



fixed points in Sec. 13.11 Consider again the Landau-Stuart oscillators 
x\ = P\X\ — uiyi + e(x\ — a) exp(^2 — (3) 

yi = Piyi + utxi. 



(10) 



where the coupling is nonlinear (and similar equations hold for subsystem 
2). The motivation for this choice of function is discussed in Sec. 13.11 and is 



common in neuronal systems [15|, [2lJ 
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Figure 8: (Color online) The largest two Lyapunov exponents as a function of coupling 
strength, at Au = (u>i = 6 and 0J2 = 6, a = 1) and j3 = 0.1, Eq. (JTDJ). Inset figure: the 
transient trajectories, x\ (solid-black line) and x-i (dashed-red line), as a function of time 
at coupling strength e = 15 [TBI l2lj. 



Shown in Fig. [8] is the largest two Lyapunov exponents for identical os- 
cillators as a function of e, the region where all Lyapunov exponents are 
negative corresponds to amplitude death (the transients are shown in the 
inset). We have found that the phenomenon of amplitude death when sys- 
tems are coupled through nonlinear interactions is very general and occurs in 
the absence of parameter mismatch, in the absence of time-delay (although 
time delay can enhance the effect) [15|, [76 . 



2. 6. Linear augmentation, and other strategies 

A final scenario we discuss dispenses with the need that the interacting 



oscillators be similar. Recently Sharma et al. [77J have proposed a new strat- 
egy through which a nonlinear oscillator, when coupled to a linear system, 
experiences amplitude death. Consider such a Landau-Stuart oscillator 

x = Px — uy + e{u — x) 
y = Py + ujx 

u = —ku — e(x — l3). (11) 

where e is the coupling strength. The variable u describes the dynamics of 
the linear system with decay parameter k (details are given in Sec. 13. 1 . 1 [) . 
and the largest two Lyapunov exponents in Fig. [9l clearly shows the region 
of amplitude death. 

Such interaction has been termed linear augmentation, in the sense that 
the dimension of the oscillator is increased by one through the addition of 
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Figure 9: (Color online) The largest two Lyapunov exponents as a function of coupling 
strength at Aui = = 6 and ll>2 = 6, K = 1) and (3 = 1, Eq. ([TT1) . Inset figure: 
transient trajectory, i, as a function of time at coupling strength e = 2 [77] . 



the linear equation. This interaction is also dynamical (cf. I2.4j) . the main 
difference between the two strategies, Eqs. (Q and ( JTTj) is that the interaction 
of each oscillator is via a linear evolution equation in the former case, while 
a single oscillator is used as self-feedback in the latter case. Such forms of 
interaction, Eq. (flT]) . are useful in the context of targeting unstable fixed 
points, as shown in Sec. 13.1.11 

Apart from above discussed scenarios, other situations where amplitude 



death occurs include the case of indirect coupling |49|, |78|, |79(: when two 
oscillators are coupled via a third, the presence of the intermediate system 
causes an effective "transmission" delay, which then effects amplitude death. 
It has been also suggested [78j that the AD is due to competition between 
synchronization and anti synchronization. Similarly, the conflict between at- 
tractive and repulsive (i.e. negative and positive) diffusive coupling also gives 
AD has also been seen in parametrically modulated systems 

and by forcing 



rise to AD 80 



8JJ , phase repulsive communication [82 



coupling [65 



83| or gradient 



2.7. Velocity Coupling 

The discussion above has focused entirely on coupled dissipative systems. 
The Lyapunov spectrum of conservative Hamiltonian systems is constrained 
by symmetry considerations to have an equal number of positive and negative 
exponents and to sum to zero. Thus these cannot show AD. However, when 
two conservative systems are coupled by time-delayed velocity coupling these 
can show AD 84 . 
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Consider, for instance two non-integrable Henon-Heiles systems 42] with 
time-delayed velocity coupling as 

oti = —x\ — 2xxy x + e(x 2 (t - r) — ±i(t)) 

y\ = -yi -xl + Vi 

x 2 = — x 2 — 2x 2 y 2 + e(ii(t — r) — x 2 (t)) 

m = -y2 + vl- (12) 

As is well-known, in the uncoupled case namely for e = the dynamics 
has both regular and irregular behavior largely depending on the total energy 



Ei although there is also considerable initial condition dependence [42j . The 
individual energies are 



Ei = + + x i y i -—. (13) 



Shown in Fig. [TOt a) are the Poincare maps for two different initial condi- 
tions, one leading to regular motion (outer points), while one leads to chaotic 
dynamics (inner points), at the same energy E = 0.13, just below the disso- 
ciation limit, E = 1/6. 

When velocity coupling is introduced say for e = 0.1 and r = 2, the 
system effectively becomes dissipative and the dynamics is attracted to the 
origin. Transient trajectories are shown in Fig. flOT b). The loss of energy 



has been discussed in a related context by Wang et. al [85| who showed 
that in the neighborhood of the fixed point, amplitude death occurs when 
the averaged total power is negative definite. In the case of the coupled 
Henon-Heiles system the energy loss of the individual oscillators, Ei and 
their difference \Ei — E 2 \ are shown in Fig. HDTc) The question of 

whether other fixed points with nonzero asymptotic energies can be targeted 
so that the amplitude death reaches a nontrivial steady state remains open. 

The above examples demonstrate that with a proper choice for the inter- 
action, a variety of systems can show AD, thus making such steady states 
amenable to control. 

2.8. Routes to AD 

In this section we consider the mechanisms for amplitude death, namely 
the "routes" followed as parameters in the system are varied. Since the 
transition to AD is from oscillatory motion to a fixed point, an immediate 



question is whether there is a bifurcation to AD 86], in an alogy w ith other 



standard routes that occur in nonlinear dynamical systems [87|, [88|, |89|, |90 . 



18 



-0.4 v 0.4 




100 1000 . 10000 



Figure 10: (Color online) (a) Poincare map at x = of uncoupled conservative Henon- 
Heiles system, Eq. (fT2]l at energy E = 0.13. Outer red points and inner black points 
correspond to the regular and chaotic motions respectively, (b) The transient trajectories 
of x\ (solid-black line) and x-i (dashed-red line) as a function of time at coupling strength 
e = 0.1 and time-delay r = 2. (c) The dissipation of energy E\ (solid-black line), E2 
(dashed-red line) of individual oscillators and their energy difference \Ei — Ei\ (dottcd- 
blue line) as a function of time. Inset figure show the repeated crossing of energies while 
decaying (Ref. 0). 
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Figure 11: (Color online) Maxima of X\ and X2 (overlapped solid lines) and real part of 
the largest eigenvalues (dashed line with circles) as a function of coupling strength, e; see 
Eq. ©. 



2.8.1. H op f bifurcation 

The most widely observed route to AD is through the Hopf bifurcation, 
where the coupling induces stability of the fixed point of the uncoupled sys- 
tems. In the Landau-Stuart oscillators, Eq. ([6]), the uncoupled system has 
two attractors, namely the stable limit cycle and an unstable fixed point. At 
the transition that occurs at e —1 (see Fig. Hj) the real part of the largest 
eigenvalue at the origin, Z* = becomes negative: this is the standard Hopf 
bifurcation since in fact a pair of eigenvalues cross the axis from right to left 



in the complex plane [88(. The amplitudes also go to zero at the same point, 
namely e = 1, as shown in Fig. [TT] for different frequencies, oji — 4 and 
u 2 = 8. 

2.8.2. Saddle-node Bifurcation 

Another route to AD when new fixed points are created is via a saddle- 
node bifurcation. These new fixed points annihilate the periodic orbit, thereby 
causing the oscillations to stop. 

For illustration consider conjugate coupled Landau-Stuart oscillators (cf. 
Eq. (ED) 

it = PiXi - ujxyx 

Vi = -Pi 2/i + wixi + ex 2 . (14) 
The coupling, which is nondiffusive results in the creation of new a fixed 
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point [111 ], in addition to the origin (0,0,0,0), at (x*, y{, x%, y^) where 



y-i 



±J(u/e(l±Veu- 



0J< 



±Jl- {x\) 2 ± Veto - u 2 , 



(15) 



with x\ 



-x\ and y\ = —y*. Note that this second fixed point depends on 
the coupling strength, and exits only when e > u (see the filled and open 
circles in Fig. [T27 a)). The trivial fixed point (0,0,0,0) (open diamond in Fig. 
IT27 b)) which exists for all values of coupling strengths is always unstable, but 
the stability of the nontrivial fixed point changes with coupling. The stable 
and unstable points collide at e = u, causing the periodic motion: thus 
amplitude deathcan be achieved with such a saddle-node bifurcation 



86 



Note also that the oscillators settle on different fixed points that depend upon 
initial conditions and thus this is an instance of inhomogeneous amp litude 
death. This route to AD has been found in other systems as well 



86 



The largest Lyapunov exponent, shown in Fig. [T2T b) shows strong fluc- 
tuations in the region prior to the AD transition, namely the marked box 
e G [cu, (1 + uj 2 )/uj} suggesting that there is multistability with both periodic 
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motion, (Ai = 0) and AD (Ai < 0) coexisting. Details are given in 
This indicates that apart from the nontrivial fixed point solutions there is 
also a periodic solution. Multistability with a riddled basin has been seen in 
diffusively coupled Rossler chaotic oscillators 11| as well as dissimilar oscilla- 
tors (Rossler and Lorenz systems) j 



SO 



In the riddled region, vanishingly 
small changes in initial conditions leaoMo different attractors, making the 
system completely unpredictable ^ 



91 



AD there is the possibility of complex dynamics [111 . I26L 1861 . 192 



This show that near the onset of 

as well. 



2.8.3. Direct transition 

Apart from the Hopf and saddle-node bifurcations, a third type of tran- 
sition to AD has been reported. Consider the system, Eq. f fT0|) where the 
interaction between the subsystems is nonlinear. Shown in Fig. [T3T a) is the 
largest Lyapunov exponent at higher value of coupling (extension of Fig. [8]): 
there is a jump from periodic motion (Ai = 0) to AD (Ai <0). This abrupt 
change in the dynamics could be a bifurcation but has not been analyzed 
completely. Similar transitions have also been observed in coupled chaotic 
oscillators as for instance two Lorenz oscillators with time-delay coupling 26 

dxi(t) 



dt 
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Figure 12: (Color online) (a) The x*i component of fixed points: trivial (blue open 
diamonds), nontrivial unstable (red open circles) and stable (black filled circles) with 
coupling strength e for conjugately one way coupled limit cycle oscillators, Eq. (JUJ), for 
uj = 2. (b) The largest Lyapunov exponent as a function of coupling strength. Marked 



box shows the multistable region where periodic as well as AD co-exist [llj, [8 
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dyijt) 

dt 
dziit) 

dt 
dx 2 (t) 

dt 

dt 
dz 2 (t) 

dt 



-x x zi -yi + riXx + e[y 2 (t - r) - y^t)] 
xiVi - pzi 
-a(x 2 - y 2 ) 

-x 2 z 2 -y 2 + r 2 x 2 + e[yi(t - r) - y 2 (t)] 
x 2 y 2 - pz 2 . 



(16) 



Shown in Fig. [TBT b) are a few of the largest Lyapunov exponents for 
a = 10, ri = r 2 = 28, p = 8/3, showing that all exponents become negative 
abruptly at a critical delay: hyperchaotic motion damps to a steady state. 



Such a transition is also observed in environmentally coupled systems (93 



but since normal form analysis is yet to be carried out, the nature of the 
bifurcation remains open. 

3. Targeting and Control 

Over the last few decades the control of both periodic and chaotic os- 
cillations in dynamical systems and the stabilization of unstable dynamics 
have been topics of intense research interest from both theoretical and ex- 
perimental po i nts o f view 87|, |94|, |95|, |96|, |97|, |98(. Several existing methods 
[9J, |95|, |96j, |99j, Il00| stabilize fixed points by changing accessible internal pa- 
rameters of a given system. In many natural systems, internal parameters are 
typically not accessible or at any rate cannot be tuned. Thus employing am- 
plitude death, using appropriate interactions between the coupled systems, 
suggests itself ClS db strategy to effect control. 

Very general methods are available for the design of specific fixed points 
in coupled nonlinear oscillators 15|, |76|, and some of these can be also applied 
to guide dynamics onto specific limit cycles. In this Section we will review 
studies that target fixed points which either exist in the uncoupled systems 
or are newly created. A related objective, the avoidance of amplitude death, 
the "anti-control" issue so to speak, is also discussed here. 

3.1. Steady-states through nonlinear coupling 

As discussed in Sec. I2.5[ nonlinear interactions can stabilize fixed points 
15l . |76| . The strategy and application of this method to achieve new and 
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Figure 13: (Color online) The largest Lyapunov exp onent as a function of coupling 
strength for coupled limit cycle oscillators, Eq. ([6]) [15(, (b) The largest three Lyapunov 
exponents as a function of time-delay for coupled chaotic oscillators Eq. (|16|) at fixed 
coupling strength e = 0.5 (iij . 
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desirable fixed points is discussed here. 

Consider coupled oscillators with a general nondiffusive coupling function 
Gj (see Eq. ([T])). The essence of the procedure to target a fixed point is the 
following. Given a set of desired fixed points, Xj, these will be stationary 
points of the coupled system with an additional constant source, namely of 
the modified dynamical system 

dXt/dt = F,(X,) + eG,(X) - F^X*). (17) 

The source function Fj(Xj) takes a constant value that depends on the 
desired fixed points. For suitable Gj it can be arranged that Gj(X) = 0. 
Upon variation of the coupling parameters such as the coupling strength e 
(or by including time-delay r in G«), the new fixed point can be stabilized: 
this, effectively, is targeted amplitude death. 

As an i llustr ation we consider coupling between identical chaotic Rossler 
oscillators [1011 ] with exponential coupling Gj = [(x, — (3) exp(xj — <5),0,0] T 



151 ] where T denotes the transpose. The resulting equations for the coupled 



system are 

ii = -yi - zi - e(xi - f3) exp(x 2 - 5) + (y x + z x ) 

yi = xi + (iy x 

zi = b + z 1 (x 1 - c), (18) 

(with corresponding equations for the other subsystem). 

The parameters and 5 in Gj are introduced such that Gj = for 
Xi = (3 and examination of the dynamical equations gives the fixed points 
Xi = (3,y~i = —(3 1 a and z^ = —bj ((3 — c),i — 1, 2. The stability of this fixed 



point can be examined as a function of e and (3 (43|. Stable (S) and unstable 
(U) regions are indicated in Fig. [TU the unstable solution corresponds to 
unbounded motion, while the stable region corresponds to the possibility of 
the AD solution. For (3 = 1, the real part of the largest eigenvalue Re(X) [43[ 
at the fixed point, shown in Fig. [15] indicates that amplitude death occurs 
when Re(A) becomes negative. Similarly the transients shown for different 
values of (3 for fixed e = 0.05 in Fig. [ToT b) arrive at different fixed points 
x~i = (3. A simple extension of this idea shows how an arbitrary point in the 



phase-space can be stabilized [76 . 
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Figure 15: (Color online) (a) Real part of the largest eigenvalue (circle with solid line) as 
a function of the coupling strength e at /? = 1. (b) Transient trajectories of x\ for j3 = 0.8 
(dashed-red line) and 1.1 (solid-black line) at coupling strength e = 0.05 [76j |. 
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3.1.1. Linear Augmentation 

A recent scheme (7?J (see Sec. 12 .6 j) to stabilize fixed points of nonlinear 
systems by coupling to a linear dynamical system, U 

X = F(X)+eU, 

U = -kU-e(X-B) (19) 

has some advantages. Since the m-dimensional linear system has the dynam- 
ics U = —kU, for positive k, and in the absence of coupling to the nonlinear 



system, this is incapable of having sustained oscillations [T7|, |78| • The addi- 
tional parameter B in the augmented system thus adaptively drives the X 
dynamics to the fixed point B. 

This scheme is ill ustra ted by the stabilization of arbitrary fixed points in 
the Lorenz oscillator 102 



x = a(y — x) + eu, 

y = rx — y — xz, 

z = xy — pz, 

u = -ku-e(x-b), (20) 

by coupling it to a linear system as above. The fixed points when e— are 
(x = ±a/ p(r — 1), = xq, z = r — 1). For a the usual set of parameter 
values, a = 10, r = 28, and p = 8/3, these fixed points are unstable, and the 
system dynamics is chaotic. 

Shown in Fig. [16] is a schematic phase diagram in the parameter space 
b — e, the shaded regime (O) representing oscillatory and chaotic motion, but 
now with an (unshaded) AD regime, corresponding to the stabilized fixed 
points. In order to achieve this we set b = x = a/ p(r — 1) = 8.4853 and 
took the decay constant k = 0.01. The largest Lyapunov exponent Ai is shown 
as a function of the coupling strength in Fig. [T7Ta) where the transition from 
chaotic motion to periodic motion and eventually to the desired fixed point 
(xq = yo = a/ p(r — 1), zq — r — 1). A typical transient trajectory of nonlinear 
oscillator in the AD regime at e = 6.8 is shown in the inset. 

Similar results are also observed for other fixed points; see Fig. [16] where 
AD occurs over a wider range of b. In the AD regime the coupling term in 
linear subsystem system vanishes since e{x — b) — > as b — )• xo, reducing the 
dynamics to u(t) = exp(-kt) namely a decay to zero. This implies that after 
a transient, the coupling term in oscillatory system will also vanish hence 
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Figure 16: (Color online) Phase space diagram in parameters (b, e) for k = 0.01 for 
b = xq. Stable fixed point solution exists in the region marked AD, while O represents 
the unstable steady state solution (oscillatory motion) ff% . 



stabilizing the unstable fixed point xq (see the inset of Fig. 1171) . Thus any 
unstable fixed point can be targeted using an appropriate value of b. 

By choosing b ^ xq, as the coupling terms do not vanish the system now 
becomes effectively 4-dimensional (in general, the augmented system has the 
dimension n + m). In this case, there are new fixed points, and since the 
region of AD in Fig. [16] is quite wide, these can also be stabilized. Of course, 



this is distinct from targeting; see details in [77 



3.2. Reprieve: Avoiding fixed points 

Some situations demand that the fixed point solution be unstable, and 
either chaotic or periodic solutions should be targeted. Although it is easy to 
perturb the system by external stimuli to avoid amplitude death, getting a 
specific solution is not always simple. Here we discuss methods that specifi- 
cally attempt to find particular periodic solutions through targeted nonlinear 
interactions; see Sec. 13.11 

A modification of the coupling function 67, for instance making it {xi — 



(3 sin(wt)} in Eq. ( fT8l) [76[ results in the dynamics being on a periodic orbit of 
frequency u. This behaviour is robust, namely it occurs in a specific range of 
parameters. If one plots the frequency of the synchronized Rossler oscillators 
as a function of the forcing as in Fig. [181 one can see that the common 
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Figure 17: The largest Lyapunov exponent as a function of coupling strength, s. Inset 
figure shows the transient trajectory, x, as a function of time at coupling strength e = 6.8, 

Eq. dm) [73. 



frequency Vt can be made to assume any desired value. Note that the other 
parameters are fixed so that in the absence of periodic forcing, amplitude 
death would result. The input and output frequencies are identical, showing 
that one can indeed target periodic motion as desired. The inset figures in 
Figs. [18] show the time series {x\ vs time) associated with such a targeted 
periodic motion for u> = 5. 

As indicated earlier, this method ensures only that the frequency of the 
targeted periodic solution takes a specific value; controlling the amplitude of 



the oscillations remains an open problem [76 



3.2.1. Gradient coupling 

The gradient (or directional) coupling introduced recently 103j | for iV >2 
is a new and interesting way to avoid AD. Note that amplitude death in 
time-delay coupled systems occurs in "islands" in parameter space (the pa- 
rameters in question being the delay time r and the coupling strength e; see 
Sec. 12.2ft . In order to avoid such islands of AD asymmetry in the coupling 
strength, keeping the coupling function the same has been used. Consider 
the modification in Eq. (J7J) as 

Zi = [1 + zcji - |Zi| 2 ]Zi + (e + r)[Z 2 (t - r) - Zt], 

Z 2 = [l + iu 2 -\Z 2 \ 2 ]Z 2 + (e-r)[Z 1 (t-r)-Z 2 ]. (21) 

The asymmetry parameter r is introduced in the coupling so that when r = 
the system is homogenous, while when r = e the coupling is anisotropic, 
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Figure 18: Variation of the common frequency of coupled Rossler oscillators (j3 
e = 0.01) as a function of targeted frequency u> [76| . 
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leading to one- way coupling [1031 ] . Such coupling arises in many practical 



103 



situations 

Shown in Fig. [19] is a schematic phase diagram in r — r space for fixed 
coupling strength e — 2. For r = 0, the system behaves as shown in Fig. 
|5] where AD is observed within the interval r G [0.248, 0.27]. When r is 
increased, the islands of AD get truncated (although the behaviour is non- 
monotonic) showing that AD can be avoided in thr ough interactions that 



have directi onali ty. A detail analysis can be found in [103 



Konishi [104] ] has shown that when the Jacobian matrix at a fixed point 
of an isolated system has an odd number of real positive eigenvalues, then 
time-delay connections are not able to induce amplitude death to this fixed 
point. It is clear that such a property can be designed, possibly through 
augmentation of the system (Sec. 12.61 & 13. 1 . 1[) into higher dimensions. This 
latter approach and the gradient method discussed above are of use when 
there is time-delay coupling. It should be emphasized that there is no gen- 
eral method that can be used in order to avoid amplitude death from any 
scenario, and thus a techniques for avoiding AD are necessary, and the prob- 
lem is worthy of further study. 



4. Characterization 

Since the asymptotic dynamics in the region of amplitude death is always 
stationarity, the interesting dynamics are essentially the transient behaviour. 
Since this is also always a result of interaction between two or more systems, 
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Figure 19: The schematic phase diagram in parameter space t — r for fixed coupling 
strength e = 2, Eq. (gTJ) . 

other phenomena that arise in coupled dynamics such as synchronization, 
riddling, and multistability can (and often do) occur, and should be taken 
into consideration in understanding the nature of dynamics near fixed points. 
In this section we focus on the nature of transient dynamics through analysis 
of the Jacobian matrix at the fixed point. AD can only occur when all 
eigenvalues have negative real parts, but the nature of the transients depend 
on whether there is a complex component or not. 

We first consider mismatched Landau-Stuart oscillators (Sec. I2.ip with- 
out time-delay. Shown in Fig. l20T a) is the frequency of individual oscillators 
as a function of the coupling strength keeping other parameters as in Fig. 

Straightforward analysis of the characteristic equation for the eigenvalues 
after setting A = a + i/3 and separating real and imaginary parts gives, for 
the leading eigenvalue, the result 



a 



l-e 




(22) 
(23) 



if e < (uji — u-2)/2, while 



a 




(24) 
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(3 = (u 1 + u 2 )/2 



(25) 



when e > (u>i — uj 2 )/2. 

This result holds when the direction of "rotation" of both oscillators is 
the same, which here corresponds to u± and u 2 having the same sign i.e. if 
one oscillator (in Fig. Ha)) are rotating in clockwise then another follows 
the same direction. When the signs differ, the anal ysis is more involved, and 



for the case Ui = —u 2 = oo (see details in Ref. (47|) the real part of the 
eigenvalue is unchanged (Eqs. fl2"3"|) and f[2"5"|) ) but 

= V(" 2 - e 2 ) (26) 

for e < oj, while (3 = if e > u. 

Shown in the insets of Figs. 1207 a) and (b) are the transients for both 
cases, namely co- and counter-rotating coupled oscillators. Fig. I2"07 c) shows 
the occurrence of AD when the signs of Ui differ; compare the behaviour 
seen in Fig. H] for the same extent of parameter mismatch. In Fig. [2"0T a) 
note that there is no synchrony for e = (cui — u 2 )/2, but for counter-rotating 
oscillators, the motion is synchronized for e < u> while beyond this value 
of the coupling the system is overdamped. These results suggest that the 
nature of the transients thus depends on the scenarios (Sec. [2]) as well as on 
whether the transition is from a globally synchronized state or not. 

4-1. The phase-flip 

When identical oscillators are coupled with time-delay (Sec. 12. 2p they 
first synchronize in phase. Shown in Fig. l2~lT a) is the common frequency and 
in (b) the phase difference between the oscillators. Within the AD region, 
there is a point at which there is an abrupt change in the frequency and this 
is accompanied by a jump in phase difference by tt: the transient trajectories 
(inset figures) go from being synchronized in-phase to being synchronized 
out-of-phase behavior. 

This transition, termed the phase-flip has been studied in detail in Refs. 



26l . 1271 . |28| . [70, [67j in a variety of systems, as well as in different dynamical 
regimes. Furthermore, when there are more than two coupled oscillators, 
the phase-difference is not necessarily ir, but can depends upon th e spe cifics 



of the system as well as on the coupling topology of the network [1051 ] . For 



asymmetric time-delays, even in the case of two coupled oscillators, the phase 



difference can be different from ir [67j. Such different phase relations have 
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Figure 20: (Color online) Frequencies, fii and fl2, for respective oscillators for (a) mis- 
matched, u)\ = 4 and U2 = 8, Eq. (|3|) and (b) counter rotation, u>± = 2 and = —2 (Sec. 
(|2.ip ~). (c) Largest three Lyapunov exponents for counter rotation. Inset figures show 
the transient trajectories in respective regimes as a function of time at coupling strengths 
e = 1.5 and e = 2.2. 
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Figure 21: (Color online) (a) The common frequency of oscillations, (b) phase difference 
between oscillators for time-delay coupling, Eq. The inset figures show the transient 
trajectories in respective regimes as a function of time at time-delays r = 0.24 and r = 
0.28. 
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been seen in experiments with asymmetric coupling Jl06[ and time-delay 67 
as well as in in relay-coupled nonlinear oscillators [79] . 

As we have seen in Sec. 12.8.21 where amplitude death coexists with pe- 
riodic motion, multistability can occur prior to the onset. The possibility 
of coexistence of both chaotic dynamics as well as stationarity has been ob- 
served in mixed Rossler and Lorenz systems 26| where there is also riddling. 



Although multistability is not a property of fixed point here, it should also 
be emphasized that at a given set of parameters reaching the fixed point 
solution can depend on initial conditions as well. 



5. Amplitude death on networks 

Most of the systems discussed in the previous sections are of two coupled 
systems only. However, many systems of interest are composed of large num- 
bers of interacting oscillators — particularly in biological contexts — which are 
coupled on networks. Mathematically a network is a graph where N nodes (or 
vertices) are connected in a specific manner. The details of the connections 
are an input in G (see Eq. ([T])). There is considerable interest in networks 
in diverse areas of the scien ces, spurred by the realization that these ideas 



have wide applicability [107]]. The internet, food webs, neural networks, elec- 



trical power grids, coauthorship an d citation netw orks of scholars, biological 



networks are well known instances 107. 108. 109 



The nature of the complex dynamical behavior on networks has a sensitive 
dependence on the nature of the connections between the nodes in addition 
to the actual dynamics of the individual subsystems. The occurrence of 



AD in networks of coupled oscillators [29|, |30|, |3l|, |32|, [33|, |34], l35| has been 



investigated for a variety of topologies ranging from small-world connections 

AD as an 
and 



38| to scale-free networks [39j and the ring topology [36|, |37 



emergent phenomenon has been seen in networks of coupled maps [55 



in chaotic coupled map lattices [61 



Since one important area for the occurrence of AD is in neuronal networks, 
we st udy here a system composed of identical Hindmarsh-Rose (HR) neurons 



llOj that are synaptically coupled in different topologies, 

e Kl 



2 3 

(XdC A $A 



iji = (a + b)x- - yi 
Z{ = cidxi + e — Zi 



(27) 
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where 

fa ~ Vs) 
[l + e X p{-/?(^(t-r)-e s )}] 



x h r) = M i „„ r K ih 77 ^ Q w (28) 



The topology is specified through the adjacency matrix elements, Aij, which 
takes the value 1 if the nodes are connected, and if they are not. At 
each node % one has a neuron, and the notation is as follows: Xi is the 
membrane potential and t/i and z\ are the fast and slow currents respectively. 
The external parameters are the synaptic coupling strength e and delay r. 
Here the reversal potential V s is fixed at V s = 2 at which the synapse is 
excitatory, the spiking threshold is fixed at Q s = —0.25 and the synaptic 
coupling function is taken to be sigmoidal. The other para meter s are set to 
standard values, a = 2.8, b = 1.6, c = 0.001, d = 9 and e = 5 [ill 



Networks of such oscillators coupled without time-delay (namely for r = 
0) continue to show global amplitude death as does the case of two oscillators. 
Fig [22] is bifurcation diagrams for networks of iV = 10 oscillators coupled in 
different topologies. Three cases considered for illustration include (a) global 
connection, namely Ki = N-l in Eq. (|2"7|) and all A^— 1 for i ^ j, (b) 
a linear chain with periodic boundary conditions, Ki = 2 as well as (c) 
random connections with Ki = 3. The motion goes from being oscillatory to 
amplitude death as the coupling strength is increased in each of the cases, 
although the threshold depends on the coupling topology. In general, though, 
the approach to the amplitude death regime is gradual with the range of 
oscillation becoming narrower with increasing coupling. 

The behavior of large ensembles of coupled weakly nonlinear oscillators 



has been studied by Ermentrout and coworkers [32j in the finite large iV as 
well as the iV — > oo limits. The state of uniform rest can become stable 
when the coupling is sufficiently strong and the frequencies are sufficiently 
disparate. Amplitude death is also reported in limit cycle oscillators co upled 
via mean field coupling and randomly distributed frequencies js), 0, [lo[ 112 



When the distribution of frequencies has sufficiently large variance, the os- 
cillators effectively pull each other off their limit cycles into the origin which 
is a stable equilibrium point for the coupled systems. 

A closed chain of delay coupled identical limit cycle oscillators that are 



near supercritical Hopf bifurcation 36( also show AD. The model is a discrete 



dynamical system, and a dispersion relation valid for an arbitrary number 
of oscillators is obtained in the continuum limit. A detailed linear stability 
investigation of these states has been carried out in order to delineate their 
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Figure 22: (Color online) Bifurcation diagrams for the system of N= 10 coupled oscillators 
with (a) global coupling, (b) nearest neighbor coupling with periodic boundary conditions, 
and (c) a small-world network where each neuron is randomly connected to two others 
(shown in the insets). x m corresponds to the maxima of x\ in Eq. (|27[) (see for details 
Ref. [H). 
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actual existence regions and to determine their parametric dependence on 
time delay. It is shown that island regions of AD can exist for any number 
of oscillators so long as the delay is finite: in effect, time delay contributes 
toward their stability. The size of these island is independent of N when N 
is even and is a decreasing function of N when N is odd. Sen and coworkers 



36l l50l . |52[ have extensively investigated amplitude death in delay-coupled 



networks. 

By applying averaging methods to diffusive time-delay coupled networks 



of weakly nonlinear oscillators [30(, it is possible to det ermi ne parameter 



ranges in which the network experiences amplitude death |ll3j . Time-delay 



also suppresses oscillations in network of oscillators coup led in a one-way 
ring topology [37|. When two networks are coupled [1 131 ] delayed effects of 
shortcuts in each network and coupling between the two groups are consid- 
ered. When the short-cut delay is removed the death of oscillations is caused 
by variational coupling delay. This behavior is found to hold in large number 
of locally coupled oscillator as well. 

As already indicated, AD in neuronal networks is a phenomenon of prac- 
tical interest, and there is some indication that such dynamics is signi fican t 
in understanding the temporal activity appearing in the olfactory bulb 114 . 
In this context there has been an analytic study of neuronal models on a 
network with time-delayed self connection. The effect of delay on the sta- 
bility of the trivial solution, an d on the existence of self-sustained periodic 
solutions has been investiga ted 114j . Neural networks with distributed delay 
also show amplitude death 115j . When the delay distribution has dispersion 
below that of the exponential distribution, the system exhibits reentrance: 
stability is lost and then recovered when the mean delay is increased. If 
the delay dispersion is greater than that of the exponential distribution, the 
system never destabilizes. 

Partial amplitude death — effectively a chim era — has been seen in a small 



network of globally coupled oscillators (116J. Konishi and coworkers [14 . 



371 . l73l . |74j have studied a number of models with different types of connec- 



tions, on networks of varying topologies. Analytical estimation of regimes of 
amplitude death have been carried out, and their analysis shows that the 
odd-number property that is known in delayed feedback control also exits in 
global dynamically coupled oscillators. Furthermore global amplitude death 
is experimentally demonstrated in a Chua's circuit coupled by an RC line. 
Various transitions including amplitude death have been also investigated in 
a regular array of nonlinear oscillators 31]. 
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A systematic procedure for achieving amplitude death involves adjust 



ing the coupling strength and delay times in the connections [61] . Such a 



procedure has the advantage that stabilization can be achieved independent 
of the network topology. This has been applied with success to the well- 
known double-scroll circuit on a small-world network (a realization can be 
seen in Fig. |22T c)). Hou et al. 



investigated oscillator death behavior on 
small-world networks and found that small-world connectivity can eliminate 
the AD seen for the same model on a regular lattice. On the other hand, 
the small-world connectivity can also lead to global oscillator death which 
is absent in the regular lattice or for that matter on a completely random 
(Erdos-Renyi) network. On scale-free networks of nonlinear oscillators AD 
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has been investigated both numerically and analytically in 

Netw orks of coupled scalar maps also show amplitude death as shown 
by Atay [l 17 ] who studied such a network with weighted connections taken 
so that they may include time-delay. The emergence of A D in a chaotic 
coupled map lattice with irregular network topologies 1181 ] or on networks 



with globally dynamical interaction 73j has been investigated in recent work. 
Similarly, amplitude death that results from conjugate coupling has been 
verified theoretically in a network of chaotic oscillators under local and global 



conditions 68 



The Belousov-Zhabotinsky reaction in a continuous-flow stirred tank re- 
actor (CSTR) is a well-defined nonlinear oscillator that can be realized in 
experiment. In addition to the study of theoretical models, experi men ts on 
networks of real systems like these also report amplitude death 119] . In 



the next Section we discuss several experimental systems where AD can be 
realized. 



6. Experiments and applications 



The stabilization of steady states, either those corresponding to the un- 
coupled oscillators or those created by design has been observed in various 
experiments. This has a variety of potential applications, and we discuss 
these here. 

Early experimental dem onst rations of AD were in coupled Belousov- 
Zhabotinskii (BZ) reactions 120| . namely chemic al os cillators. Other exper- 



iments on asymmetrically coupled BZ reactions |106] also lead to AD. The 



experimental set up used identical BZ reactors coupled via mass transfer- 
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the coupling is thus diffusive in nature and since it depends on the concen- 
tration gradient both symmetric and asymmetric coupling can be realized by 
changing the active volumes of the tank reactors. This is of importance in 
biological systems where systems are indeed coupled through the exchange 
of fluid and are likely to be asymmetric in nature. A situation that is partic- 
ularly amenable to experimental control is that of electrochemical oscillators 
121] ; by coupling nearly identical systems in the proximity of their individ- 



ual Hopf bifurcation points, they can be driven to AD |6j as discussed in Sec. 



Thermo-optical oscillators (122J, where the coupling is linear and by heat 



transfer are ano ther class of coupled systems that show AD. Thermo kinetic 



oscillators [1231 ] that show temperature and concentration oscillations are 
of direct relevance in industrial processes. Most of these are exothermic in 
nature, and AD, which has been seen experimentally in such systems can 
have severe conseq uence s; these studies thus may have a significant impact 
on safety measures 
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Many systems and scenarios for AD have been experimentally investi- 
gated in coupled electronic circuits. Apart from the Landau-Stuart system, 
the van der Pol oscillator and other standard nonlinear or chaotic oscillators 
that are easily realized in analog circuits, a number of hybrid models can also 
be studied in experiment. This extends the study of mismatched oscillators 
beyond the cases of dy namic coupling [37J or when the frequencies of coupled 



units differ [124j, Il25j to the situation when the coupled systems are com- 
pletely distinct from each other: the interaction between an electronic and 
biological oscillator — a ne ural-electronic hybrid interface — can also suppress 



oscillatory dynamics (126 



Partial amplitude death is the situation when inst ead of all variables 
reaching a stationary state, only some of them do 1271 ] . This behaviour is 



seen when strong interactions break down both reflection and translational 
symmetry in the coupled syst ems, and has been experimentally verified using 



coupled Lorenz circuits [128|. Partial AD is distinct from the case of the 



chimeric behaviour in networks, where some members of a set of coupled 



osci l lator s remain oscillating while others go to a stationary state [3l|, [55 



129, 130 



Experiments on coupled lasers have been a major area for the applica- 
tion of many ideas discu ssed above since delay coupling can be easily imple- 
mented. Pyragas et al. 13l| have proposed that the stabilization of fluctua- 



tions in multi-mode intracavity doubled diode pumped Nd:YAG lasers can be 
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achieved by stabilizing the steady state, namely through AD. A r elated ex 



perimental study of chaos control in lasers by Bielawski et al. [1321 ] stabilizes 



the unstable steady state of a fiber two-level class B laser. Delayed cou pling 



with delayed feedback has been shown to induce AD in a laser system [133 



and Roy and coworkers |l6l . 1 17] have experimentally verified a number of 



issues in the dynamics of coupled systems, including amplitude death. 

The fact that several meteorological phenomena are oscillatory and are 
frequently coupled suggests that AD will also occur in climate models. There 
are large-scale oceanic and atmospheric decadal anomalies. This variability 
is characterized by patterns of sea-surface t emp erature anomalies, correlated 
with the atmospheric pressure at sea level 134] . In the ocean, this variabil- 
ity results in fluctuations of ocean temperature and ocean mass transport. 
In the atmosphere these fluctuations are associated with changes in surface 
westerlies and surface air temperature over ocean basins and the period of 
these oscillations is appears to depend on the transit time o f oce anic Rossby 
waves; these are of the order of decades. Gallego and Cessi 134] . using their 
model for mid-latitude large scale interaction between the upper ocean and 
the troposphere consider the case when two ocean basins are coupled through 
zonally average atmosphere. Each ocean basin also interacts with the atmo- 
sphere via wind driven torques and heat fluxes at the sea surface. When 
uncoupled, each individual ocean basin has sustained oscillations, but with 
coupling, the oscillatory anomalies decay in time, and the system eventually 
reaches a steady state. 

This vital result in the study of climate change is thus linked to amplitude 
death which occurs only when the long Rossby wave delays for the two ocean 
basins are sufficiently different. The model formulated here uses the Rossby 
wave delay times to control the natural frequency of oscillation in each basin, 
and thus the case of AD seen here can be seen as an example of mismatched 
interacting oscillators as discussed in Sec. 2.1. 

Other large-scale systems that show AD include ecological models as for 
example coupled prey predator systems 135] . Amplitude death is important 
in epidemical models as well, e.g. in infectious disease, its important to pre- 
dict whether the infection has disappeared or the pathogen persists. Steady 
state s in the models of pathogen-immune dynamics of infectious diseases 



136] are other instances of systems where the AD state is a desirable end 
point, and thus methods of control, discussed in Section 4 are important in 
the context of intervention strategies. 

The issue of amplitude death has been also discussed in the field of neural 
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oscillators (See Eq. (I27p a nd corresponding Fig. |2"2"|) . Amplitude death in a 
neural network model |l37l | where the interaction is "transferable" , namely in 
a cyclic chain of oscillators with one-way coupling and with time lag. Ampli- 
tude death has also been observed in a network of Fitzhugh-Nagumo neurons 



138]. Diffusive coupling is not the only means of achieving amplitude death. 



Phase-repulsive coupling [82J also eliminate oscillations in a population of 
synthetic genetic clocks. 



7. Summary 

The phenomenon of amplitude death has attracted considerable attention 
in the past decades following the realization that this is quite general and 
that this has significant practical implications. 

In this review we have shown that AD is a general outcome in coupled 
nonlinear oscillator systems, occurring in a wide variety of dynamical models 
and for different forms of coupling. There are a number of scenarios which 
lead to amplitude death and we have tried to use the well-studied Landau- 
Stuart limit cycle oscillator model to illustrate these. A detailed theory — in 
terms of bifurcation structure — is not yet available, and this constitutes an 
interesting and open problem. 

The asymptotics of AD is a featureless steady state, but the transient 
dynamics can be of interest, with the motion of the coupled systems having 
nontrivial phase relations. The frequency of the damped oscillations can also 
show a discontinuity: this is the phase-flip transition. We have reviewed 
here these various dynamical characteristics in model systems. 

Targeting a specific steady state is a central issue in the study of nonlinear 
dynamics, particularly with regard to stabilizing low-dimensional dynamics 
such as fixed points or periodic orbits. We have reviewed some aspects of 
targeting fixed point solutions by considering specific interactions. This will 
be essential if one is to be able to use the phenomenon of AD in practical 
applications. Targeted fixed power outputs in coupled laser systems can have 
significant applications in laser technology e.g. laser surgery, or laser welding 
and fabrication. In other practical applications AD is to be avoided, for 
instance in brain function; we have also reviewed methods for the reliable 
avoidance of AD. There can be significant constraints, though since many 
natural systems have specific and constrained coupling forms — for instance 
the synaptic coupling between neurons, or the diffusive interaction between 
nonlinear chemical systems. 
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Although we have considered simple systems for purposes of illustration , 



AD is known to occur in higher dimensional systems as well |139l . Il40l . Il41 



Time-delay makes the dynamics effectively infinite-dimensional and thus 
there can be AD in hyperchaotic systems. 

External noise is unavoidable in natural systems or in experiments. AD 
is quite robust with regard to noise since the Lyapunov exponents of systems 
undergoing AD are substantially negative; the systems are deep in the region 
of stability. Thus AD has been observed in numerous experiments many of 
which have been discussed here. This aspect deserved further study, however, 
since it is not clear to what extent inhomogeneous steady states can be 
perturbed by noise (e.g. after e = 2.5 in Fig. IT2"]) . In addition, the study 
of network is currently of great importance in various areas of study and as 
we have seen, AD in complex systems with different topologies is of interest. 
The chimeric state of partial amplitude death is, in particular, an important 
one that needs further exploration. 

Although there has been substantial work on amplitude death, several 
challenges remain both at fundamental as well as applications levels. Through 
this review we have tried to summarize the current state of understanding of 
the phenomenon, as well as drawn attention to several open issues, and to 
highlight possible applications. 
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